Testing Gravity using the Environmental Dependence of Dark Matter Halos 
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In this Letter, we investigate the environmental dependence of dark matter halos in theories which 
attempt to explain the accelerated expansion of the Universe by modifying general relativity (GR). 
Using high-resolution iV-body simulations in f(R) gravity models which recover GR in dense envi- 
ronments by virtue of the chameleon mechanism, we find a significant difference, which depends on 
the environments, between the lensing and dynamical masses of dark matter halos. This environ- 
mental dependence of the halo properties can be used as a smoking gun to test GR observationally. 
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One of the biggest challenges in cosmology is to ex- 
plain the recently observed accelerated expansion of the 
universe. The acceleration might originate from either 
"dark energy" within the framework of GR, or from a 
large-scale modification to GR without introducing new 
matter species. It could be difficult to distinguish between 
these two scenarios by merely measuring the expansion 
rate of the Universe, and one has to study the growth 
of structure formation in the Universe to break the de- 
generacy. On large scales, it is possible to perform model 
independent tests of GR by combining various cosmolog- 
ical observations [1-3], but information on linear scales 
is limited due to theoretical degeneracies as well as sta- 
tistical and systematic uncertainties in observations. 

There is ample information available about cluster- 
scale structure formation, but it is difficult to predict the 
observablcs on non-linear scales in modified gravity (MG) 
models. If GR is modified on large scales, there may ap- 
pear in gravity new scalar degree of freedoms, (dubbed 
scalaron), modifying GR even on cluster scales. In order 
to evade the stringent constraints on deviations from GR 
in the solar system, we need a mechanism to recover GR 
on small scales by screening this scalar mode. In such a 
mechanism, e.g. the chameleon mechanism in f(R) grav- 
ity [4, 5], the mass of this scalar mode depends on the 
local density of matter, becoming heavier in denser envi- 
ronments and thereby suppressing the scalar interaction. 
This environmental dependence can provide us with a 
smoking gun for alternative theories to GR [6, 7]. 

In this Letter we investigate and quantify for the first 
time the environmental dependence of the difference be- 
tween lensing and dynamical masses for dark matter ha- 
los. Our analysis is based on our high-resolution iV-body 
simulations for the f(R) gravity model [8], where the 
Einstein-Hilbcrt action in GR is extended to be a gen- 
eral function of the Ricci scalar [9] . 

In the Newtonian gauge in a general perturbed 
Friedmann-Roberston- Walker universe, the line element 
can be written as ds 2 = a 2 (?7)[(l + 2$)d?7 2 - (l-2*)df 2 ] 
where r\ is the conformal time, a(r\) is the scale factor, 



$ and VP are the gravitational potential and the spatial 
curvature perturbation respectively. The Poisson equa- 
tion reads 



V 2 $ 



47rGa 2 <5peff, 



(1) 



where G is Newton's constant, and Sp c g is the perturbed 
total effective energy density, which contains contribu- 
tions from matter and modifications to the Einstein ten- 
sor due to MG. The dynamical mass Mrj(r) of a halo is 
defined as the mass contained within a radius r, inferred 
from the gravitational potential felt by a test particle at 
r. It is given by Mr> = J* a 2 Sp c ffdV , in which the integral 
is over the extension of the body. Under the assumption 
of spherical symmetry, the Poisson equation can be inte- 
grated once to give 



M D (r) oc r 2 d$(r)/dr. 



(2) 



To measure Mo from our N-body simulation we use the 
force acting on particles to infer the force acting on each 
halo function of the halo radius and Mr> can be 

obtained using Eq. (2). Observationally, Md can be esti- 
mated from measurements such as velocity dispersions of 
galaxies. In f(R) gravity Md includes the contribution 
from the scalaron, which mediates the finite-ranged fifth 
force within the Compton wavelength. The mass of the 
scalaron depends on the local density of matter, resulting 
in the environmentally-dependent modifications to Md- 
On the other hand, the lensing mass is determined by 
the lensing potential $+ = ($ + \&)/2. In f(R) gravity 
for example, <3>+ satisfies V 2 <E>+ = AirGa 2 8pM, where 5pM 
is the matter density fluctuation if we assume that the 
background cosmology is close to that for ACDM. This 
is the same equation as in GR, since the scalar mode 
does not couple to photons and it does not modify light 
propagation [12]. The lensing mass is defined as Ml = 
J a 2 SpMdV, and is the actual measured halo mass in our 
simulations. Thus we will use Ml to represent the halo 
mass throughout. For a spherically symmetric body we 
have 



M L (r) oc r 2 d$+(r)/dr. 



(3) 



2 




9 9 



9 Halo mass M , log 1Q D>l 
9 Halo mass M L , log lQ D<l isolated halos, log lQ D>l clustered halos, log l0 D<l 



llog A I llog A I 

I &io M 1 1 & 10 M 1 



FIG. 1. The illustration of the three-dimensional distributions of the halos (panel A) and the mass difference Am for each halo 
(B,C) for the \f R0 \ = 10~ 6 model. The sizes of the bubbles are proportional to the halo mass Ml in (A), and |log 10 Ajvf| in 
(B,C). In all the panels, the pink and blue bubbles illustrate the isolated halos, and the halos living in the dense environment 
respectively. See Eq. (5) for the definition of D, which quantifies the environment. 



The lensing mass and the dynamical mass are the same 
in GR, but they can be significantly different in MG sce- 
narios. To quantify the difference, we calculate the rela- 
tive difference A m between Ml and Mo for each halo, 
Aj\/ = Md/Ml — 1. Similar quantity, g = Am + 1, was 
introduced in Ref. [11]. Combining Eqs (2) and (3), we 
can rewrite Am as, 



. , . d$(r)/dr 



(4) 



In GR Ajvf(r) = 0, while in MG models Am{t) varies 
depending on the local density. 

We have chosen f(R) gravity as a working example 
to investigate how Ai/(r) correlates with both the halo 
mass and the environment, and propose a new method 
to test GR based on this correlation. For the analysis 
we shall use the high-resolution iV-body simulation cata- 
logue [8] for a f(R) gravity model, f(R) = aR/(/3R + j) 
[10] where a = —m ci,0 = C2,7 = — m 2 ,m 2 = HqCIm 
and ci,C2 are free parameters. The expansion rate of 
the universe in this f(R) model is determined by C1/C2, 
and the structure formation depends on |/ro|j which is 
the value of \df/dR\ at z = 0, and is proportional to 
ci/c|. We tune c\/c2 to obtain the same expansion his- 
tory as that in a ACDM model, and choose values for 
|/flo| so that those models cannot be ruled out by cur- 
rent solar system tests. To satisfy these requirements, 
we set c\/c2 = QQ^/SIm and simulate three models with 
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FIG. 2. The contour plots between halo mass Ml and D for 
f(R) and GR models on a log-log scale. The shaded colour 
stands for the number density rescaled by the average number 
density of halos in each pixel on the Ml-D plane. 
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SR(fn) where 



perturbation of the Ricci scalar 6R(fn) = 
_ 3V^<5/h anc j - g fluctuation of fa = 
can see that when the scalar mode vanishes, 
0, we recover the GR relation between cur- 
matter, SR = SitGSpm and Am = 0. This 
the dense region where the chameleon has 
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FIG. 3. AM(r34o) as a function of halo mass (panels 711-^44) and D (B1-B4) for our f(R) models on a log-log scale. Each 
magenta dot represents one halo in each case, and the black solid line and the error bars show the mean value of Am (r 340) and 
the l-<7 error bars respectively. The correlation coefficient p is shown in all cases. 



an effect, but in the underdense region where 6R can 
be ignored, 5p e g = |<5pm so Am = 1/3. One expects a 
strong correlation between Am and Ml since halos with 
large Ml should be 'screened' against the modified grav- 
ity influence, and GR be locally restored. This has been 
confirmed by the previous analysis [8, 11]. 

The mass threshold for the screening can be estimated 
theoretically [11]. Interestingly, we find that the small ha- 
los with masses below the screening mass threshold can 
also be well screened if they live in dense environments. 
This effect is shown visually in Fig. 1. In panel (A) we 
show the 3-D map of the halo distribution in the f(R) 
model with \fn.o\ = 10~ 6 where the size of the bubbles 
is proportional to Ml, while in panels (B, C) the bubble 
size is proportional to | log 10 Am | ■ In other words, a larger 
bubble means a more massive halo in (A) , while it means 
a better screened halo, in which GR is better restored, in 
(B, C). In all the panels, the pink and blue bubbles illus- 
trate the isolated halos (logio-D > 1), and the halos living 
in the dense environment (logio-D < 1) respectively, and 
these two subsets of halos are complimentary (See Eq. (5) 
for the definition of -D, which quantifies the environment). 
Panels A and B look almost the same in pattern, mean- 
ing that more massive halos are better screened, and the 
halo mass is the only factor affecting the screening. This 
is natural since the environmental effect is removed in 
panel B by design. On the other hand, in panel C, halos 
living in dense environments are shown, and the environ- 
mental effect is the dominating factor for the screening, so 
that halos with mass below the screening mass threshold 



can also be efficiently screened. The difference between 
the panels (A, C) indicates a clear environmental depen- 
dence of Am - small halos can be well screened by their 
neighbouring halos. 

This implies that Am correlates with not only Ml, 
but also the environment. The environment effect was 
also noticed by Schmidt [11]. In this Letter, we shall 
quantify this effect for the first time. The environmen- 
tal dependence of Am provides valuable information for 
testing GR, which compliments the information of the 
mass-dependence of Am ■ The amount of information can 
be maximised if the estimates of the halo mass and envi- 
ronment are uncorrelated. 

The 'environment' can be defined such that it suits the 
physical set-up of the problem, facilitates ease of obser- 
vations, or both [13]. For our purpose, we need an envi- 
ronment indicator which can represent the local density 
well, but with least correlation with Ml- Such a quantity 
was found in Rcf. [13], 

L>NJ = , (5) 

?'NB 

which is defined for a halo with mass Ml as the distance 
d to the Ath nearest neighbouring halo whose mass is 
at least / times as large as that of the halo under con- 
sideration, rescaled by the virialised radius tnb of that 
neighbouring halo. Clearly, a large value of D^j indi- 
cates a scarcity of nearby halos, meaning that the consid- 
ered halo lives in a low-density environment. It is found 
that in GR, D\ \ is almost uncorrelated with the halo 
mass, and represents the local density well [13]. 



4 



To test the mass- independence of Di i in the context 
of modified gravity, we select the resolved halos from our 
high-resolution f(R) and GR simulations with boxsize 
B = 64 Mpc/h [8]. In our simulations, the halo mass is 
measured using Ml = 4ir x Np cr ar^/3 where rjv is the 
radius when the density reaches N times of the critical 
density of the Universe p C rit, and we choose N = 340 
[8]. To be conservative, we only select the well- resolved 
halos from our simulations, i.e. halos more massive than 
10 12 /i _1 M©. In Fig. 2, we show the contour plots between 
Ml and D, (we will use D to represent D\ i hereafter 
for brevity), for three f(R) models in comparison with 
the ACDM model simulated using the same initial con- 
ditions. The darkness of the shaded colour quantifies the 
number density of the halos in each pixel on the M^-D 
plane. We follow Ref. [13] to use the Spearman's rank 
correlation coefficient p, which is the correlation coeffi- 
cient between the ranked variables and varies from —1 to 
1, to quantify the correlation between Ml and D. As we 
can see, they are almost uncorrelated in all cases since 
the absolute value of the correlation coefficient p is much 
less than unity. This means that the information of the 
A]\i-D relation is highly complimentary to that of the 
Am~Ml relation, which provides us with a new means of 
testing gravity observationally. 

Fig. 3 shows Aa/ (Y340) as functions of Ml and D for 
two f(R) models. To further disentangle the residual cor- 
relation between Ml and D, we divide the samples into 
three subsamples in both cases. In the A panels the ha- 
los are divided according to their ordered D values. Halos 
with D values in the top third of the group (logio-D > 1) 
are classified as halos in an 'Underdense Environment' 
(Al,A3), while those with D values in the lowest third 
(logio-D G [—1,0.65]) are viewed as halos in an 'Over- 
dense Environment' (A2,AA). In the B panels, the halos 
are separated according to their mass, namely, the halos 
whose mass is in the top third (Ml > 10 12 ' 7 Mq//i ) are 
called 'Large Halos' (B2,B4), and the third with small- 
est mass (M L G [10 12 , lO 12 - 3 ]M /ft.) are labeled as 'Small 
Halos' (Bl,B3). The horizontal blue dashed line shows 
AA/(f34o) = 1/3, which is the threshold of Am in f(R) 
gravity. 

As can be seen from Fig. 3, for the |/i?o| = 10~ 6 model, 
Am(^34o) decreases when Ml increases, as expected. 
Note that this anti-correlation is stronger (p = —0.74) in 
the underdense regions (panel Al), as the environmental 
effect can be safely ignored in these cases, and Am(?"34o) 
is mainly determined by Ml- In the overdense environ- 
ment (A2 for the \fno\ = 10 -6 case) for the halos with D 
values in the lowest third, the effect of external environ- 
ment becomes important - many halos less massive than 
10 12 ' 3 Mq//i get screened thanks to the boundary condi- 
tions set by neighbouring halos. An interesting observa- 
tion is that some small halos are better screened than the 
big ones in this case. This is because many small halos 
reside in overdense environments, surrounded by many 
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FIG. 4. The profile of logioAju as a function of the rescaled 
halo radius r/r 340 for the |/ho| = 10 -6 model. We show the 
profile with 1 — a error bars for the halos divided into four 
categories as illustrated in the legend. The red dashed line 
shows A M (r34o) = 1/3. 



neighbouring halos or inside very big halos, while large 
halos are more likely to be isolated so that their screening 
is mainly determined by their mass. For the |/ho| — 10~ 4 
model, the halos are very weakly screened in all cases, and 
Am(^34o) is close to 1/3, which is the maximum relative 
mass difference in f(R) gravity. The \fno\ = 10~ 5 case is 
somewhere in between, and is not shown here. 

The environmental effect can be seen more clearly in 
the A M (r 340 )-D plot (panels B1-B4 in Fig. 3). For the 
I /no I = 10~ 6 case, we see a strong correlation between 
the two, and this correlation is largely insensitive to Ml 
(p ~ 0.7) in both mass bins. Again, we see that the 
screening is very efficient in dense regions even for the 
least massive halos. The correlation between Am(?"34o) 
and D for the \fm\ = 10~ 4 case is much weaker, which 
is because essentially none of the halos are screened by 
either their own masses or the environment. 

We show the profile of the mass difference as a func- 
tion of rescaled radius for the |/ro| = 10~ 6 model in Fig 
4. To see the environment effect on the profile, we split 
the samples according to both the halo mass and D pa- 
rameter. As we can see, small halos in the underdense 
region are hardly screened at all, while the halos with 
similar mass in the dense region are efficiently screened, 
and the screening effect is stronger in the core of the ha- 
los. For large halos, the innermost part is well screened 
regardless of external environment due to the high mat- 
ter density there, but the part close to the edge shows a 
clear environmental dependence, and the difference can 
be as large as 3 orders of magnitude in Am in different 
environments. This is because in this region the external 
environment plays an important role. 

The lensing mass and the dynamical mass can be mea- 
sured using strong lensing and the peculiar velocity dis- 
persion measurements, respectively, and there has been 
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some effort to test GR by comparing the two observation- 
ally [14, 15]. However, the measurements of the absolute 
values of Am are likely to be contaminated by system- 
atics. Fortunately, the strong environmental dependence 
of Am due to the scalar mode in modified gravity theo- 
ries may provide a way to ameliorate this problem. Ob- 
servationally, one could divide the galaxy samples into 
different groups using D, and measure the difference of 
Am among those subsamplcs. If a Am~D correlation is 
found, then it can be viewed as a smoking gun of a mod- 
ified gravity signal, which can be independently tested 
using the Am-Ml correlation. 

In this Letter, we focus on the Chameleon mechanism 
to recover GR on small scales. There are different classes 
of mechanism to achieve the screening, such as the Vain- 
shtcin mechanism [16] and the symmetron mechanism 
[17]. In the case of the Vainshtein mechanism, it was 
found that the screening of halos is almost independent 
of the environment [11]. Thus the method we proposed 
provides not only a new independent test of GR on fully 
nonlinear scales but also a way to distinguish between dif- 
ferent screening mechanisms. It is extremely interesting 
to perform this test using the high-quality observational 
data from the upcoming large-scale structure surveys. 
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